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We calculate the distribution of the conductance G in a one- dimensional disordered wire at finite 
temperature T and bias voltage V in a independent-electron picture and assuming full coherent 
transport. At high enough temperature and bias voltage, where several resonances of the system 
contribute to the conductance, the distribution P{G{T, V)) can be represented with good accuracy 
by autoconvolutions of the distribution of the conductance at zero temperature and zero bias voltage. 
The number of convolutions depends on T and V . In the regime of very low T and V , where only one 
resonance is relevant to G{T, V), the conductance distribution is analyzed by a resonant tunneling 
conductance model. Strong effects of finite T and V on the conductance distribution are observed 
and well described by our theoretical analysis, as we verify by performing a number of numerical 
simulations of a one-dimensional disordered wire at different temperatures, voltages, and lengths of 
the wire. Analytical estimates for the first moments of P{G{T,V)) at high temperature and bias 
voltage are also provided. 



PACS numbers: 72.10.-d, 73.23.-b, 72.15.Rn 

I. INTRODUCTION 

The continuous progress in the fabrication of small 
electronic circuits has kept active the theoretical study 
of quantum electronic transport. For example, transport 
properties in point contacts, atomic chains, carbon nan- 
otubes, and quantum wires are currently under exper- 
imental investigations^iSiSi^i^. Although the theoreti- 
cal study of electronic transport in mesoscopic systems- 
where the phase coherence of electrons is preserved along 
the whole device- has been of great interest for several 
decades, these recent experimental advances have further 
renewed the motivation of theoreticians for studying the 
electronic transport through one-dimensional structures. 

In a disordered mesoscopic wire the random position of 
impurities gives a stochastic character to the electronic 
transport. Therefore, it is of particular relevance the 
analysis of the statistical properties of transport quanti- 
ties like the conductance. On the theoretical side, there 
is a full description of the statistical properties of the 
conductance in quasi-one-dimcnsional disordered systems 
within an independent-electron picture at zero tempera- 
ture T and "zero bias" (actually, infinitesimally small 
voltage) Vi^ This degree of detail in the theoretical de- 
scription is, however, not available for disordered wires at 



finite temperatures and bias voltage, even in the simple 
case of one-dimensional (ID) wires. This is an unfortu- 
nate fact since experiments are usually performed within 
a wide range of T and V . For example, the effect of finite 
bias voltage has been found to modify the behavior of the 
conductance fluctuations^^ 

Several efforts have been made in order to incorporate 
the information of a finite T into the statistical descrip- 
tion of quantum electronic transport in ID-disordered 
systems. In Ref. jfl, at zero voltage, a model of zero- width 
resonances represented by (5-functions was introduced in 
order to calculate the dependence on the temperature of 
the averages of the conductance {G{T)) and (In G'(T)). 
Assuming full coherent transport, in Ref. JJ3 the effect 
of the thermal smearing on the mean resistance of a ID 
wire was studied. 

At finite temperature and bias voltage, the distribution 
of the conductance P(G(T, V)) was calculated in Ref. O 
by using a statistical model of resonant tunneling trans- 
mission. The methodology there employed was, however, 
restricted to disordered wires of length L ^ I, where I 
is the mean free path, and very small values of T and V 
(although sizable): typically eV and kT were considered 
of the order or smaller than the mean spacing between 
energy levels A. Within this regime of temperatures and 
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voltages it can be assumed that only one resonance-the 
closest to the Fermi energy- contributes to the trans- 
port. Even with these limitations, that work shows that 
the conductance distribution do display novel features as 
a consequence of finite T and V values. Particularly, it 
was shown that the distribution of conductances narrows 
as the value of T and/or V is increased. 

We would like to remark that the assumption of phase 
coherent transport adopted in previous works, as well 
as in the present one, is restricted by the phase relax- 
ation at finite temperature and voltage. For exam- 
ple, inelastic scattering and thermally activated hopping 
processeSfi2iiii4 which might be important at finite val- 
ues of T, y are not considered in our analysis. In Ref. 
ITU however, it is discussed the possibility of satisfying 
> L even at finite temperature and bias voltage for 
short enough wires, which opens the possibility of observ- 
ing the effect of finite T, V on the statistical properties of 
the conductance in realistic systems. It is also interesting 
to note that the range of voltage V where the averaged 
charge current {J{T, V)) is a linear function of V or de- 
parts slightly from this behavior, i.e., the so-called linear 
response regime, can be rather wide as we will show in 
Section Irnl 

The main goal of the present paper is to compute the 
distribution of the conductance G'(T, V) at values of eV 
and kT typically larger than the mean level spacing, 
where several resonances contribute to the conductance. 
We show that, within this regime of temperatures and 
voltages, the distribution of the conductance can be well 
described in terms of the convolutions of the zero temper- 
ature/voltage conductance distribution p{g). The num- 
ber of distributions to be convolved is determined by the 
number of resonances in an energy window where electron 
transport can take place. This energy window is defined 
by the value of the temperature and voltage. The method 
introduced in the present work applies to any degree of 
disorder i.e. to any value of the ratio L/l, provided sev- 
eral resonances contribute to the transport, allowing for 
the investigation of short wires where coherent transport 
is more likely to be observed. For completeness in the 
sense that all regimes of temperature and voltage will be 
covered in this paper, we also study the regime of small 
T and V following Ref. ,11 

This paper is organized as follows. In Sec. ^ we 
present the methodology to analytically calculate the 
conductance distribution P{G{T^V)). We divide this 
section into two subsections. The first one, III Al is 
devoted to review the case of small temperatures and 
voltage regimes where one can assume that only the 
closest resonance to the Fermi energy contribute to the 
conductanceii In the second one, subsection lll Bl we in- 
troduce a method based on the convolution of the known 
distribution of conductances at zero temperature and 
bias voltage in order to study the regime of high tem- 
peratures and bias voltages. 

In section IIIII we present the model used in our sim- 
ulations and show some general features of our ID- 



disordered system obtained numerically, which support 
the assumptions introduced in our theoretical proposal 
of section m In section Hill we also compare the results 
of the distribution P{G{T, V)) from the numerical simu- 
lations to the predictions of the theoretical approach pre- 
sented in the previous section, for the different regimes 
of T and V. Finally, in section IIVI we give a summary 
and conclusions of our study of the distribution of con- 
ductances at finite temperature and bias voltage. 

II. METHODOLOGY 

We consider the usual setup where the disordered con- 
ductor is placed between left and right reservoirs at 
the same finite temperature T and chemical potentials 
Mi = M + eV/2 and fiR — fi ~ eV/2, respectively. The 
conductance G{T,V) = J{T,V)/eV can be written as 
(in units of the conductance quantum 2e^//i) 

G(T, V) fdEgiE) (f{E - eV/2, T) -f{E + eV/2, T)), 

(1) 

where f{E, T) = {exp[(£' - ^Ji)/kT] + 1} Ms the Fermi 
function, being k the Boltzmann constant, while g{E) is 
the dimensionless conductance for E = fj, a,t T = 0. It is 
clear from the above integral expression for G{T, V) that 
the energy window where electronic transport takes place 
will depend on the values of T and V. This can be em- 
phasized recasting the difference in the Fermi functions 
in the above expression, Eq.([3l, asi^ 

fiE-eV/2,T)^ fiE + eV/2,T) = 

{Q{E-eV/2)-Q{E + eV/2)) * ^M|ll),(2) 

where the symbol * denotes the convolution in energy, 
<d{x) is the unit step function, and df{E,T)/dE is the 
thermal smearing function. Therefore, the difference of 
the Fermi functions can be expressed as a convolution 
in energy of two functions: one depends only on the ap- 
plied voltage V and the other one, only on temperature 
T. This implies that thermal and voltages effects are sta- 
tistically independent, and this will be useful to calculate 
the distribution of G(T, V). 

A. Small temperatures and bias voltages: 
one-resonance contribution to G{T, V) 

In this subsection we briefly review the resonant model 
introduced in Ref. fTll This model is appropriate for 
systems with non overlapping resonances, which is sat- 
isfied for disordered systems with L > I. Assuming a 
Lorentzian line shape for the resonances, the dimension- 
less conductance g{E) at T = 0, y ^ can be written 
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where = tI > + is the total width, tI '''' cx 
Aexp[— (L±2zi,)/^;^] are the left and right partial widths, 
and are the localization lengths. On resonance {E = 
E,y), the terms of the sum Eq. — [cosh(2z^/^y)]~^, 

depend on the location z,y of the state, and are maximum 
for z,y = (center of wire). Substituting Eq. |j2J into Eq. 
PJ, it is found that the conductance G(T, V) is given by 



1. Zero temperature, finite bias voltage 

Let us first consider the simple case of finite bias volt- 
age and zero temperature. In this case, the thermal 
broadening function in Eq. (j^J is a delta function and 
therefore g{E) in Eq. is only multiplied by a rectan- 
gular function of width eV. Thus Eq. (jSJ is reduced to 
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where Wn = (2n + l)7rfcT. We now assume that << 
kT, eV^A. Under this condition only the resonance clos- 
est to the Fermi level {v — 0) contributes to the conduc- 
tance G{T, V) and therefore Eq. Q can be simplified 
toil 



tanh 



Eq- ^1- eV/2 
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2kT 



(5) 



In order to calculate the distribution of the conduc- 
tance P{G{T,V)) a statistical model for the resonances 
is introduced: the resonances Eq and their positions zq 
are considered uniformly distributed in a interval A and 
L, respectively, while the distribution of the inverse of 
the localization lengths p(a;o — 2L/£_o) follows a Gaussian 
distribution with mean (xq) ~ L/l and var(a;o) — 2{xq). 
This distribution for p(xo) is good for systems with 
I <^ L. In Section IlIII we will show some examples for 
the distribution P{G{T,V)) obtained within this reso- 
nant tunneling model. 



B. Large temperatures and bias voltages: 
contribution from several resonances to G{T,V) 

We now go to the main goal of this work and propose 
a description of the conductance distribution, which is 
valid for arbitrary length of the wire and for voltages and 
temperatures satisfying eV, kT > A, but small enough to 
satisfy that the behavior of the average current does not 
depart from linear response. 

Firstly, we discretize the integral in Eq. ^ as 
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-f{E,-ifi + eV/2),T)], 



(6) 



where we have assumed the same width SE for all the 
elements in the sum, Eq. © . 



(7) 



i=l 



where the number of terms in the sum A'' satisfies NSE= 
eV. This approximation is exact in the limit N oo. 
However, being our aim the evaluation of P{G{T,V)), 
we approximate by a finite number N of statisti- 
cally independent contributions. We now assume that 
this number corresponds to the mean number of levels in 
an energy window eV, i.e., N = eV j IS.. This is a natural 
assumption since for a given disorder realization g{E) is 
a spectral function with peaks at the energy levels of the 
wire (resonances). In a non- interacting electron picture, 
the height of the resonance peaks, as well as the energy 
levels change for different disorder realizations, but the 
average separation between levels is a well defined quan- 
tity and the spectral weights centered at the different en- 
ergy levels are uncorrelated. In addition, we assume that 
g(E) is a random stationary function of the energy, at 
least in the energy window where transport takes place. 
Then the statistical properties of g{E) do not change 
in such energy window and therefore the distributions 
p{g{Ei)) are actually independent of the energy Ei, i.e. 
pigiE,))=p{g). 

Under the above assumptions the distribution 
P{G{T = 0, y)) can be computed from the convolution 
of N distributions p{g), i.e., the iVth autoconvolution of 
the distribution at zero temperature and bias voltage: 



P{G{T 



= 0,y)) =p(l)(g)*p 



(2) 



(.9)*...*pW(g), (8) 



where we have defined g = g/N , and the upper indices 
enumerate the number of distributions p{g) that enters 
into the convolution. 

As we mentioned in the Introduction, at zero temper- 
ature and infinitesimal small bias voltage the statistical 
properties of the transport quantities, in particular the 
distribution p((7) for ID and quasi-lD disordered systems 
is well known. In a framework of random-matrix theory, 
a diffusion equation known as Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equation has been successful in describ- 
ing the evolution of the conductance distribution p{g) as 
a function of the system length L in quasi-one dimen- 
sional systems^. For strictly ID wires the DMPK equa- 
tion is reduced to the Melnikov equation whose solution 
p{g) can be written asS 
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FIG. 1: (Color online) Sketch of the behavior of the function 
llljl l and its approximation by a rectangular function, in black 
dashed and red solid lines. The behavior of the integral of the 
function 1131 1 is also indicated in dashed-dotted lines. 



where i/q = arccosh(2/g— 1) and s = L/l is the length 
of the system L in units of the mean free path I. In the 
limit of s >> 1 and s << 1 closed analytical expressions 
for p{g) can be obtained&ii We remark that the only 
parameter that enters in the distribution p{g) is the so- 
called disorder parameter s = L/l. 

The simple relationship between the random variables 
G{T — 0,V) and g, Eq. {T)), encloses important con- 
clusions since this means that the Cumulant Generating 
function for G and g are also simple related: 



logAfG(A) - A^log(Af§(A)) 



(10) 



with Af^(A) = / exp(— A7) P{j)d'y being ^ — G,g, while 
for the average and the first three central moments Eg = 
((G - (G))«) and = {{g - {g))'>) it is satisfied: 
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for g = 2, 3 



(11) 
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(Ti = by definition. Then, Eq.l^TJ shows that the mean 
value of the conductance (G) is independent of the bias 
voltage V. 



2. Finite temperature and bias voltage 

Let us now consider a more realistic situation where 
temperature and bias voltage are both finite. 

At finite temperature, the difference of Fermi functions 
in Eq. |^ reads 



fiE~ip~eV/2),T) 
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In fig. ^it is shown the behavior of Eq. (|13|l . In order to 
keep our statistical analysis of the conductance as simple 



as possible, we approximate the bell-shaped function H13() 
by a rectangular function of height of tanh(ey/4fcT) and 
width eV / tanh{eV / AkT) . The area of the rectangle is, 
therefore, eV as the area of the original bell-shaped func- 
tion Eq. H13I) . This simplification allows us to proceed as 
in subsection 111 B l( from the width eV / tanh{eV / 4kT) 
we define an "effective number of resonances", N^ff, 
given by 



eV 



eff 



A tanh(ey/4fcT) ' 



(14) 



We can verify that for T ^ 0, N^ff N ^ eV/A, as 
expected. Thus to calculate P{G{T,V)) we can follow 
exactly the same procedure of the previous subsection 
111 B II we just replace N by Nf,ff. Also, at finite T and 
V it is possible to define a simple relation between cu- 
mulants and moments of the distribution P{G{T, V)) and 
p{g); again, substituting N by A^e// in Eq- l|T^ . 



III. NUMERICAL RESULTS 

In this section we verify our theoretical study of the 
statistical properties of the conductance, in particular 
the conductance distribution, by comparing to numeri- 
cal simulations of a ID-disordered wire. The results pre- 
sented here give numerical evidence that supports the 
main hypothesis of our theoretical model and benchmark 
the quality of the approximation of the distribution func- 
tion at finite bias and temperature on the basis of con- 
volutions of the distribution in Eq. ^ against exact 
numerical results. 

In our simulations we model a ID-disordered wire of 
length L using the 

standard tight binding Hamiltonian of spinless elec- 
trons with a single atomic orbital per lattice site and 
nearest neighbors hopping parameter t: 



(13) 



Cj + l 



H.c.) + SicXci 



(15) 



with L = A^sfl, being a the lattice constant, which we set 
as the unit of length and A^s , the number of sites of the ID 
lattice. The on site energies Ej are chosen randomly from 
an uniform distribution of width W (Anderson model). 
All the energy scales are taken dimensionless in units of 
the hopping parameter t. The disordered conductor is 
connected to the left and to the right to ID clean semi- 
infinite leads, which we assume to be at chemical poten- 
tials — ^J. + eV/2 and = fi — eV/2, respectively, and 
temperature T. The Hamiltonian describing the contacts 
between the reservoirs and the disordered wire reads: 



(16) 



being k^^n the contact sites of the left and right semi- 
infinite leads, respectively. 



For each disorder realization, the current is numerically 
calculated from the expression with 

giE)^r^iE)\G^^SE)\', (17) 

being r{E) ^ \t\^e{\E\ - U)y/m^ - E^ / At^ , which cor- 
responds to reservoirs modeled by semi-infinite tight- 
binding chains with hopping i, while Gfj^{E) is the re- 
tarded Green's function for the disordered chain con- 
nected to the reservoirs Typically, we generate numer- 
ically 10'* to 10^ realizations of disorder. For several val- 
ues of the disorder parameter s — L/l, we have verified 
that for very small bias {eV << A) and T = we are 
able to reproduce the distribution p{g) given by Eq. 
The mean free path I is extracted from the numerical 
simulations through the relation (Ing) = —L/l. In all 
numerical simulations we present below we have fixed 
jjL — and the disorder strength to W — 0.5 which sets 
I = 144. 

The theoretical distributions P{G{T, V)), as described 
by Eq. ||SJ) of the previous section, are obtained by col- 
lecting the data for G{T, V) from an ensemble of conduc- 
tances generated numerically using a Metropolis Monte 
Carlo algorithm. 




FIG. 2: (Color online) {J{T, V)) as a function of V ior W = 
0.5, ^ = and different lengths L = 120 (solid ), L = 200 
(dotted), L = 400 (dashed), and L = 800 (dashed-dotted). 
Inset: {G(r, V)) vs L for T = and V = 0.01 (black crosses), 
V — 0.1 (blue squares) and V = 0.2 (red diamonds). Note 
that the different symbols are overlapped. 



A. General features 

Before presenting our results for the probability distri- 
bution at finite temperature and bias voltage, we would 
like to establish the range of the bias voltage V where 
the averaged current {J{T = 0, V)) varies linearly with 
V. To this end we have included in Eq. ((TKl) the poten- 
tial drop due to the presence of the bias by modifying 
the local energies as Ej Sj + eV/2 — jeV/Ng. In Fig. 
121 we show the disorder averaged current {J{T,V)) as a 
function of the bias voltage, at T = 0. A linear behavior 
of {J{V)) with the bias V is observed up to < 0.5 for 
different lengths of the system. Therefore if we want to 
restrict our study to the linear response regime, we can- 
not take arbitrary large values of T and V. In the present 
model, this implies energy values for eV and A:T^0.5. Let 
us note, that this range of energy is actually rather wide: 
it is larger than 10% of the total band width of the clean 
system (equal to At). Since we focus in voltages within 
the linear response regime, in what follows we consider 
the wire model (Eq. (|15l) ) without including the effect of 
the potential drop in the local energies Sj. 

In the inset of Fig. [2 we also show the average of 
the conductance {G{T, V)) for different values of T, V, 
as a function of the system length L. We observe that 
{G{T, V)) is independent of the values of the temperature 
and voltage and decreases exponentially with L as it is 
expected. This behavior is in agreement with previous 
resultsii as well as with the prediction based on Eq. ((TT|l . 

Another important point to verify is the reliabil- 
ity of our assumption that the mean energy spac- 
ing between levels of Ea. H15|l sets the energy scale 



where the spectral weights of g{E) are statistically 
independent. To this end, we have investigated 
the behavior of the auto-correlation function C(e) = 
{g{E)g{E + e)) - {g{E)){g{E + e)) , where the over line 
means energy average. In Fig. 13 we have plotted C[e) 
(normalized by dividing by C(0)) for three different wire 
lengths: L = 40, 120 and 800. In all the cases C{e) de- 
cays between 90% ^ 99% of its value at e = for e < A. 
In order words, C{e) is a vanishing function in an energy 
scale e smaller than the mean level spacing. 

We have also verified numerically that in the linear 
response regime p{g{E)) is in fact independent of i?, as 
we assumed in the previous section. 

On the other hand, it is instructive to show the evo- 
lution of the conductance distribution P{G{T,V)) with 
the number of autoconvolutions of p{g) , as described in 
the previous section. We have chosen the simple case of 
T = with and finite V = 0.1 with s = 400/144 p^i 2.7. 
In Fig. ^ the histogram in dashed line corresponds to 
the distribution obtained from our numerical simulations, 
while the histograms in solid line correspond to P{G) 
obtained by m autoconvolutions of p{g), for to = 4, 6, 8, 
and 10. We can observe that as m increases from 4 to 
6 the theoretical P{G) evolves to the numerical distri- 
bution. When m = N = 8, which corresponds to the 
mean number of levels N ior L — 400, see Fig. O a very 
good agreement is found with the numerical distribution. 
The fact that the optimum number of autoconvolutions 
of p(g) coincides at T = with the mean number of lev- 
els within eV, supports the validity of our approach of 
subsection III B II 
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800 



FIG. 3: (Color online) Mean number of levels n of the Hamil- 
tonian 11511 in an energy interval equal to 0.1, as a function of 
the system length L. Inset: The autocorrelation function C(e) 
is plotted for L = 40, (A = 0.108), L = 120, (A = 0.046) 
and L = 800, (A = 0.007) (red open circles, black solid cir- 
cles and blue square symbols, respectively). See the text for 
details. 
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FIG. 4: (Color online) Histograms obtained from Eq.© per- 
forming different number of convolutions m, for a finite volt- 
age y = 0.1 and L — 400. a) m = 4, b) m = 6, c) m = 8 and 
d) m = 10. For comparison the distribution P{G) ( histogram 
in red dashed line) obtained numerically is also plotted. For 
m = 8 (= A'^) convolutions the distribution P(G) is quite 
well reproduced. The size bins of the numerical and theo- 



B. Small temperatures and bias voltage: P{G{T,V)) 
from one resonance contribution to G(T, V) 

In this subsection we show some numerical results for 
the distribution P{G) when the values of the temper- 
ature and bias voltage are small enough that only one 
resonance contributes to the conductance G{T, V) (sub- 
section III A|) ■ We also indicate the limitations of this 
method by showing an example of the conductance dis- 
tribution at regimes of T and V beyond the scope of the 
resonant model. 

The theoretical distribution P{G{T, V)) is calculated 
by generating numerically an ensemble of conductances 
G accordingly to Eq. |SJ , where the random variables Eq 
and zq are obtained numerically from a uniform distribu- 
tion, while xq is extracted from a Gaussian distribution, 
as it is described in lll Al 

In figure we show the distribution of In G (we have 
chosen the variable In G since the details of the distribu- 
tion can be better seen in the logarithm scale) from the 
resonant tunneling model (histogram in solid line) for the 
case of L = 400 ~ 31 and eV = A. A good agreement 
with the numerical simulation (histogram in dashed line) 
is seen. Note that, although this situation corresponds to 
a small bias voltage, the exact distribution (solid lines) 
clearly departs from the behavior of the "zero" bias dis- 
tribution p{g) (histogram in dotted lines). Therefore the 
resonant model reasonably describes the numerical be- 
havior of P{G), in particular the decrease of the width of 
the distribution as V,T increases (with V,T < A). How- 
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FIG. 5; (Color online) Histograms obtained with the reso- 
nant model of subsection III Al (solid black lines) for L — 400. 
Results in dashed red lines correspond to the numerical sim- 
ulations on the disordered tight-binding chain. Plots in the 
main panel correspond to V = A,T = while the inset cor- 
responds to V = A,T = A/3. The distribution p{g) is also 
shown in the main panel (blue dotted line). 



ever, if we increase the value of T, keeping eV — A, in 
our previous numerical example, the quality of the ap- 
proximation provided by this model deteriorates, as it is 
illustrated in the inset of the figure [S| This happens in 
general when the values T and/or V are increased in such 
a way that more than one resonance might contribute to 
the conductance. In the next subsection we will see that 
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FIG. 6: (Color online) Distributions P{G) for V = 0.2 and 
T = (red dashed line histogram) obtained numerically for 
different lengths a) L = 40, b) L = f20, c) L = 200, and 
d) L = 400. The solid black line histograms correspond to 
the distribution of G obtained as a convolution of A'^ terms a) 
N = 2, h)N = 4, c)iV = 7, d)iV = 14. The blue dotted curves 
show p{g) for zero bias and temperature. The size of the 
bins of the numerical and theoretical histograms are slightly 
different to distinguish better between both histograms. 



TABLE I: moments S, {q = 2, 3) for bias vohage V = 0.2, 
W = 0.5 for different lengths L obtained numerically and ct, 
from the distribution p{g), Eq.Q. 



L 


iV 


S2 

xlO"^ 


(72 /iV 

xlO^^ 


Eg 

xlO"^ 


(73 /iV^ 

xl0~^ 


40 


2 


1.469 


1.427 


-1.32 


-1.01 


80 


3 


1.716 


1.802 


-0.574 


-0.460 


120 


4 


1.554 


1.733 


0.010 


0.054 


200 


7 


1.099 


0.933 


0.32 


0.263 


400 


14 


0.414 


0.383 


0.160 


0.110 


800 


29 


0.065 


0.071 


0.015 


0.012 



TABLE IL Moments (g = 2, 3) for bias voltage V = 0.01, 
W = 0.5 and several temperatures T for L — 400 obtained 
numerically and aq from the distribution p{g), Eq.Q . 



r 


Neff 


Ea 
xl0~^ 


a2/Neff 

xlO"^ 


E3 

xlO"^ 


xlO"^ 


0.005 


2 


1.964 


2.585 


3.24 


5.0 


0.01 


3 


1.175 


1.723 


1.239 


2.2 


0.05 


15 


0.294 


0.344 


0.107 


0.088 


0.1 


29 


0.144 


0.178 


0.023 


0.023 



2. Finite temperature and bias voltage 



the convolution method describe correctly the regime of 
T and V when several resonances are involved in the 
transport problem. 



C. Large temperature and bias voltage: P{G{T,V)) 
from several resonances contributions to G{T, V) 

1. Zero temperature and finite bias voltage 

Let us now go to the analysis of the probability distri- 
bution function at finite V with T = using the method- 
ology introduced in subsection lllB II In figureElwe show 
the distributions P(G(T = Q,V)) for different lengths L 
of the disordered wire with V = 0.2. We recall that 
I = 144 in all cases. For each length L, the distribution 
P{G) obtained from Eq.© (histogram in solid line) and 
the numerical distribution (histogram in dashed line) are 
both displayed for their comparison. A very good agree- 
ment can be seen. In order to provide evidence of the 
strong effect of the finite bias voltage on the conduc- 
tance distribution, p{g) is also shown in dotted lines in 
the same figure. 

We have also verified the relation between the second 
and third moments, Eq. H12|l. For a bias voltage V — 
0.2, in Tabled we show the values of extracted from 
the numerical simulation and those for aq {T ^ V ^ Q) 
calculated from then integral expression of p{g), Eq. 
Again, a good agreement between numerics and theory 
is obtained. 



Finally, let us consider both finite temperature and 
bias voltage. This case corresponds to the regime con- 
sidered in subsection lll B 21 There, an effective number of 
resonances was introduced, Eq. p4|l . which corresponds 
to the number of autoconvolutions of the distribution 
p{g), Eq. ® , to be used in order to obtain P{G{T, V)). 
In Fig. 13 we present the results for P{G{T,V)) (his- 
tograms in solid line) from the convolution method for 
L = 400, V = 0.01, and four different values of T. For 
each value of T, the number of convolutions N^f f changes 
according to Eq. (|14|) . As in previous cases, we compare 
to the numerical simulations (histograms in dashed line). 
A good agreement is seen in all cases. 

The first moments of the conductance distribution 
were also obtained. The results are shown in Table ^1 
where, as before, is obtained from the numerical simu- 
lation, while tTq is computed from the expression for p{g), 
Eq. The agreement between the numerical simula- 
tion and theory is, in general, reasonably good. We recall 
that at finite T we have introduced the simplification of 
representing the difference of the Fermi functions by a 
rectangular function, Eq. (|13|l . This additional approxi- 
mation might be the cause for the discrepancies between 
and <Jq/NlJj observed for small values of A^e// (first 
two lines of table II). Notice that a small error in the esti- 
mate of the integer number N^f f implies a large relative 
error for small N^ff- In all cases, however, the behavior 
of the moments and the conductance distribution is well 
described by our method. 

Finally, it worth mentioning that we have verified that 
the approach based in convolutions of p{g) also provides 
a good approximation to the exact P{G{T,V)) at finite 
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FIG. 7: (Color online) Numerical distributions P{G) (dashed 
red line histograms) of a system with length L = 400, V = 
0.01, and different temperatures: a)T = 0.005, b)r = 0.01, 
c)r = 0.05, d)T = 0.1. The theoretical distribution from 
N^j f autoconvolutions of p{g) (solid line histograms) are com- 
pared with the corresponding numerical results: a,)Neff = 2, 
h)Neff — 3, cjNeff = 15, d)Neff = 29. A good agreement is 
observed. 



T in cases with larger eV (e.g. eV = 0.1, 0.2). 



IV. SUMMARY AND CONCLUSIONS 



tunneling model, Eq. |(SJ, as we have verified numeri- 
cally. As the temperature and bias voltage is increased 
several resonances contribute to G{T, V) . In this regime, 
the conductance distribution P{G{T, V)) can be obtained 
from the convolutions of the known distribution of con- 
ductance at zero temperature and bias voltage ^(17). The 
number of autoconvolutions of p(g) is determined by the 
width of the energy window where transport can take 
place. In the case of zero temperature and finite bias 
voltage, the width of the energy window is trivially eV 
and the number of convolutions is given by the mean 
number of levels N — eV/A. In the case of finite tem- 
perature and a finite bias voltage, we have simplified the 
problem by introducing an effective number of resonances 
that allow us to reduce the problem of finite T and V to 
the simpler case of zero temperature, and finite V. The 
results of our theoretical method have been compared 
to numerical simulations of a ID-disordered system at 
different regimes of temperature, voltage, and different 
values of the length of the system. A good agreement 
has been found in all cases. We point out that for small 
T and V the line shape of the resonances was relevant 
in the calculation of G{T,V) (see Section FlI All . When 
the values of T and V are such that several resonances 
contribute to the conductance, the line shape is seen to 
be irrelevant. To conclude, we remark that with the res- 
onant model and the convolution method of sections III Al 
and lll Bl respectively, we are able to analyze the conduc- 
tance distribution at all regimes of temperature and bias 
voltage, under the assumptions described in the paper. 



In an independent electron picture and assuming full 
phase-coherent electronic transport, we have studied 
the distribution of the conductance P{G{T, V)) in ID- 
disordered systems at finite temperature and bias volt- 
age. 

We have observed a strong effect of finite T and V 
at the level of the conductance distribution. In general 
P{G{T, V)) is narrower compared to the conductance dis- 
tribution at T = 0, y ~ 0. The average of the conduc- 
tance is, however, independent of T and V (see Fig. 
This implies that higher moments have to be analyzed to 
see the effect of the temperature and voltage. 

When temperatures and voltages are small (less than 
the mean-level spacing) , only one resonance is relevant to 
the conductance. In this regime, the distribution of the 
conductance is well described by a simplified resonant- 
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